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Abstract 

A family of threshold parameters which probe the stability of chiral predictions is considered. The 
relevant criteria for the choice of threshold parameters are discussed. Sum rules for these quantities 
are derived from dispersion relations and evaluated from effective range formulae. Good agreement 
with two-loop chiral estimates for many of these quantities is found and interesting discrepancies are 
discussed. 
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1. Dispersion relations for vrvr scattering amplitudes with two subtractions have been rigorously es- 
tablished in axiomatic field theory It is convenient to consider dispersion relations for s— channel 
amplitudes of definite iso-spin, T/(s,t, n), 7 = 0, 1,2, where s, t, u are the Mandelstam variables. 
Furthermore, each of the amplitudes may be written down in terms of a unique function as T^{s, t, u) = 
3^(s, t, u) + A{t, u, s) + A{u, s, t), T^{s, t, u) = A{t, u, s) - A{u, s, t), T^{s, t, u) = A{t, u, s) + 
A{u, s,t). Unitarity, analyticity and crossing symmetry have been used extensively to study this funda- 
mental process of elementary particle physics. Introducing a partial wave expansion for these amplitudes 
via T/(s, t, u) = 327ri;(2/ + l)//(s)P;((t — n)/(s — 4)), elastic unitarity implies above threshold and be- 
low the four-pion threshold that the partial wave amplitudes may be described in terms of the phase shifts 
5l{s) by //(s) = ^Js/{s — 4) exp(i5/(s)) sin^f (s), where we have set the pion mass (771,^ = 139.6 
MeV) to unity. Note the threshold expansion Re//(z^) = i^'(af + bfv + cfv'^ + dju^ + ...), v = 
(s — 4)/4 > being the square of the three momentum in the centre of mass frame, which defines 
the scattering lengths aj, the effective ranges bj and the higher threshold parameters cj and dj, etc. 
Sum rules have been established for (combinations of) scattering lengths and effective ranges in the past 
employing analyticity and crossing symmetry constraints. One culmination of the dispersion relation 
approach to vrvr scattering has been the Roy equations [Q, ^ which is a system of coupled integral equa- 
tions for partial wave amplitudes which further trades the two unknown t-dependent functions in fixed-t 
dispersion relations for the scattering lengths aj] and Aq. Roy equation fits to phase shift data [0, |5|] have 
been extensively studied. Best fits to data give Oq = 0.26 ± 0.05 [^. 

2. Chiral perturbation theory is the effective low energy theory of the standard model and describes 
processes involving pionic degrees of freedom, viewed as the approximate Goldstone bosons of sponta- 
neously broken axial-vector symmetry of massless QCD. At leading order the vrvr scattering amplitude 
is given by the Weinberg result from PCAC: A{s, t, u) = {s — l)/F^ [^. (In chiral perturbation theory 
this amplitude has now been computed to one-loop and even two-loop order |T^.) Furthermore 
this implies that the only non- vanishing threshold parameters are aj] = 7/ (327rF^), Cg = — l/(167rF^), 
al = l/(247rF^), 6j{ = l/(47rF^) and 6g = -l/{8TrF^). It is important to note the well known 
result that the set of functions: 4° = (2/3(s - 4/3) + 5a/3)/{F^), = f5{t - u)/{F^), = 
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(— /3(s — 4/3) + 2a/3)/(F^), where a and f3 are arbitrary real constants, verifies dispersion relations 
with two subtractions and vanishing absorptive parts. The Weinberg amplitude is a special case of this 
general linear amplitude, with q = /? = 1. It may be noted that a generalized version of chiral perturba- 
tion theory that is motivated by considerations of a small quark condensate in the QCD vacuum allows 
a to vary over a range between unity and as much as three, and reorganizes the chiral power counting 
[pl|]. In this discussion we confine ourselves to the more predictive standard chiral perturbation theory 
with a = 1 for much of our discussion. 

Dispersion relations with two subtractions have been used to write down sum rules for (combinations) 
of some of these threshold parameters, for example the Wanders sum rules which maybe written down 
as: 

1 P ^ d 
47r^ Jo ./UAV+iW 



dv 



4^^ Jo TR^^+iF 

6bl = -2aO + 5al + — [2a\v) + iva\v) + {7v + - ^{2v + 1)^7^(0)] , 

Jo v(^(^ + l)) 

where a^{u) = {At: / y/u{u + 1)) ■ ^(2/ + l)Im //(z^) are the cross-sections. The Weinberg predictions 
for the quantities involved in these sum rules obey these relations identically with all the cross-sections 
set to zero, since the original amplitude obeys dispersion relations with vanishing absorptive parts. Fur- 
thermore, it may be noted that since these are relations at leading order in the chiral expansion, it would 
be fair to expect that the relevant dispersion integrals make numerically less significant contributions to 
18a}, 36o and Gb^ compared to 2aQ — Soq, even at higher orders in the chiral expansion. The presence 
of two subtractions in the dispersion relations where aj] and Qq play effectively the role of subtraction 
constants, render them difficult to pin down on the basis of dispersion relation phenomenology alone: 
vital chiral inputs are required to make sharp predictions for these quantities. 

3. At one-loop order, chiral perturbation theory requires the introduction of several coupling constants 
in C^^\ which account for the non-renormalizable character of the non-linear sigma model Lagrangian, 
which is the basis of the Weinberg result. Four of these constants Ij, i = 1,2,3,4 enter tttt 
scattering and in particular Ii and I2 enter predictions for cj], Cg, b\,a2 and a\, while all higher threshold 



parameters will receive no contributions from the trees generated by C^^'^ at this order. In the past 
experimentally known values for a2 and 02 were used to fix these quantities [Q], while no values have 
been reported in the literature for c[ and b\. For completeness we note the one-loop formulas for c[ since 
they have been reported earlier: 

Recently, rapidly converging sum-rules were written down [12] in order to estimate h\ which in 
principle could be used to fix the values of the two chiral coupling constants of interest, indeed, as one 
could from the values of Cg should they be known. However, systematic ambiguities inherent to one- 
loop predictions for these quantities have also been discussed recently. Ii and I2 were fixed instead by 
rewriting the chiral and axiomatic representations of the scattering amplitudes to the appropriate order 



in the momentum expansion Q13|]. This method is now being extended to the two-loop case in order to 
fix the coupling constants that enter at that order in CHPT. Note that the trees generated by C^^^ will 
now contribute to the threshold parameters dQ,c\, h2 and a\. At the one-loop level they are given by 
the expressions: = -1643/(403207r3F^), = -893/(403207r3F^), c\ = -23/(134407r3F^), 
ftO = -481/(2016007r3F^), hi = -277/(2016007r3F^) and a\ = ll/(940807r3F^). 
To summarize, we note that the trees generated by C^'^\ contribute to 



while those generated by C^^'^ contribute to 



al (2) 



"0 "0 '^0 

a\ b\ (3) 
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and those of C^^^ contribute to 



"0 "0 ''0 "o 



a\ h\ c\ 



(4) 



4. The main purpose of this letter is precisely to provide estimates to those quantities appearing in eq.^ 
and eq.(Q) for which no information is available in the Uterature and to compare whenever possible with 
the predictions of CHPT. Such a consistency check may be viewed as a probe into the range of validity 
in energy of chiral predictions. Indeed, when Roy equation fits to the planned precision experiments are 
performed, all the quantities discussed here may be evaluated afresh, which would then amount to a high 
precision experimental determination of these numbers. Such Roy equation fits may then be employed to 
evaluate chiral parameters determined from dispersion relation phenomenology, and such a consistency 
check may then be performed again. 



5. We work in the approximation that the absorptive parts are modeled entirely by the S- and P- waves 
(In the numerical analysis we further assume that the contribution to the dispersion integrals from the S- 
and P- waves also above the KK threshold may be neglected as in Ref. [|l^].) which has been found 
to be justified phenomenologically in the past and is supported today by chiral power counting. In this 
approximation it is particularly convenient to represent the iso-spin amplitudes as ^ [I^ : 



T'{s,t,u) = 327r ^ (hsl'^' + tC'J' + nCii')ai' + 



I'=0 



dx 



1 

vr Ji x{x — 4) 
+3 



+ 



X — t 



+ 



X — u 



Im/o^'(x) 



(5) 



s(t - u)l''' ^ t{s-u)Cll' ^ u{t - s)Cii' 



Imfl'ix) 



X — S X — t X — u 

Projecting these amplitudes onto partial waves will yield the Roy equations for each of the waves in the 
S- and P- wave approximation. [This approximation is equivalent to setting the so-called driving terms of 



the Roy equations to zero. However, in the numerical work described below, we evaluate the contribution 



of the /2(1270) resonance [ ]15| ] to the threshold parameters of interest as a measure of the contribution 
of the driving terms.] The most convenient manner in which the sum rules of interest may be computed 
is to consider the Roy equations for each of the partial waves /q , f^, and /g in the neighborhood of 
the threshold, as power series in v. Furthermore, in order to isolate the quantities of interest we need to 
consider the Cauchy Principal Value of the relevant integrals. The Principal Value singularity occurs due 
to self-coupling of the waves and must be removed during the process of computing the power series of 
the real parts of the waves. An example is worked out for the / = S- wave in Appendix A and the other 
waves may be treated analogously. The complete set of sum rules for the quantities of interest that have 
not been published in the past is listed in Appendix B. 

6. The spirit of this work will closely follow that of Ref . |jT^ wherein a modified effective range formula 
was employed in order to model the absorptive parts of the amplitudes. Note that the dispersion integrals 
we encounter for (most of) the quantities of interest are very rapidly converging which implies that it 
is the near threshold region that needs to be modeled accurately, the region where the effective range 
formula is applicable. The modified effective range formula for the S- and P- wave phase shifts is of the 
type first proposed by Schenk [ p^ and is: 



tan 5^(1/) = U^[ai + [h^ - ai/4 + {aif]A 



,1 



i^n5\{v) = {a\ + [b\ - a\/ul]A (6) 



7. Our numerical work requires inputs to the effective range formulae. Much of these inputs are guided by 
one-loop chiral perturbation theory. We also perform a sample computation with a set of input parameters 
with Oq corresponding to the best fits to the experimental data of K^^ of 0.26 and the rest of the quantities 
computed for that optimal Roy equation fit. In Table 1 the complete set of inputs is tabulated, besides the 
conventional choice Vq = 8.5, vl = 6.6 , = —5.0. 

In Table 2 we present the computed values of the threshold parameters of interest for the inputs of 
Table 1 . The contribution of the /2 resonance [[l^] in the narrow width approximation to the quantities 



of interest, calculated by plugging in the appropriate absorptive part into the relevant Roy equation and 
computed in the appropriate limit, yields the results tabulated along /2 in Table 2. This should be con- 
sidered as setting the scale of the corrections arising from all the higher waves and the higher energy tail 
that would be described in terms of Regge parameterization of the absorptive parts and possible Pomeron 
contributions. It may be judged from this that the bulk of the contribution is received from the S- and P- 
wave low energy absorptive parts. 

In Table 3 we tabulate the one-loop and two-loop predictions for these quantities whenever available 
[^, 17]. Let us first consider the one loop predictions for Cg which receive contributions from li and l2- 
The latter computed from optimal Roy equation fits of Ref. [13] have been inserted into eq.(|T]) to produce 
the entries in the first two columns of Table 3. The rest of the quantities are of the pure one-loop variety 
and their values numerically tabulated. Alternatively, one may take the values of the Cg and solve for Ii 2 
via. eq.([T]). Taking the extreme values [including the contribution of the /2] for —0.0049 < Cg < 0.010 
and 0.012 < < 0.015 we obtain -4.7 <h< 0.08 and 3.8 < I2 < 5.7. While these numbers are 
not to be taken literally since the contributions of the high energy tail are completely neglected and the 
effective range formula is not a real substitute for a Roy equation solution of the lowest partial waves, 
they continue to provide an important consistency check on the values of these coupling constants. These 
are consistent with several prior determinations for these quantities. 

The comparison of the values of Cg with the two-loop predictions continue to be encouraging as 
we observe from Tables 2 and 3. Some attention may be paid to the quantity c\: at one-loop order it 
is negative where as the two-loop prediction is positive and larger in magnitude. This value creeps up 
towards the value provided by the sum rules. It could be that the presence of the p as a non-perturbative 
feature of hadron dynamics is responsible for the mismatch between even the two-loop prediction and 
the sum rule result for c\. 



Furthermore, the Taylor series of Re//(z^) = y sin(2(5/^(z^))/2 computed with the effective 
range formulae for the phase shifts, itself yields coefficients cj and dg which are comparable with the 
quantities computed from the sum-rules. This proves to be a check on the effective range parameter- 
ization itself. For instance, for the choice I in Table 1, we find Cg = 0.0092, Cg = 0.011, c\ = 
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0.00092, = -0.013, = -0.0012, whereas for that of choice IV, we have eg = -0.0025, eg = 
0.011, c\ = 0.00084, d[| = -0.015, = -0.0013. 

Thus we see that detailed considerations of unitarity, analyticity and crossing in the near threshold 
region is in excellent agreement with chiral predictions — these, however, do not suffice to discriminate 
between the standard and generalized scenarios of chiral perturbation theory. 

8. In order to come to grips with these numbers, we have also chosen to compare the predictions for the 
P- wave threshold parameters with numbers arising from resonance saturation with the p. In particular, 
the narrow width formula for the absorptive part generated by the p: 



gives rise to an amplitude A{s, t, u) by inverting eq.(||) via A{s, t, u) = (T'^(s, t, u) — T'^{s, t, u))/3, 

AnsJ.u)= W'^; + (8) 

'rripK'nip — 4)'^/^ — t m-^ — u 

This is identical to the result in eq. (C.9) of Ref. ^ obtained from effective Lagrangian techniques with 
the appropriate identification of the relevant coupUng constant when we take into account the identities, 
l/{ml -t) = 1/ml ■ (1 + t/{ml - t)) and [t{s - u) + u{s - t)] = [-2{s - 2f + l/2{s'^ + {t - uf)]. 
This is a consequence of two subtractions in the dispersion relations used here. The formula eq. (|8|) 
when inserted into the 1 = 1 amplitude and projected onto the P- wave yields for a\ - (2a!] - 5ao)/18 = 
0.0073, b\ = 0.0077, c\ = 0.0006, when parameters of the resonance [ p^ ] are inserted into the formulas 
for these quantities. Note the p contribution to li and I2 are —7rmp/{2Tp) and 7Tmp/{ATp) respectively 
[see, e.g.,eq. (C.IO) in 



Another treatment of resonance saturation | |18[ ] with the p is based on writing an unsubtracted disper- 
sion relation for the 1=1 t-channel amplitude divided by (s — n). This is expected to converge under the 
assumption of the validity of the Pomeranchuk theorem, supported by the behavior of the relevant Regge 
trajectory. The p contribution to li and I2 are, however, —TTmp/{4rp) and 7rmp/{4Tp) respectively [as 
read off from eq. (28) and eq. (32) of [jlSl]]. This is in variance with eq. (C.IO) of [^]. This must result 
from considering dispersion relations without a sufficient number of subtractions that guarantee the ab- 
sence of crossing constraints on the absorptive parts due to the P- wave. [An analogous situation arises 
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with the resonance saturation with states of Z > 2 of fixed-t dispersion relations with two subtractions; 
this does not yield a crossing symmetric amplitude and must be treated with care.] 

9. In the preceding sections we have estabhshed a program of comparison between families of threshold 
parameters for which sum rules have been derived and the quantities numerically estimated, and then 
compared with one- and two-loop predictions for these quantities. The role of the number of subtrac- 
tions in dispersion relations is crucial and has been emphasized in a variety of ways. Other recent efforts 
in this confrontation and comparison between data, dispersion relations and chiral perturbation theory 
use extrapolation to subthreshold regions reaching analogous conclusions regarding the interplay of an- 
alyticity and unitarity [|^. Our work is also in keeping with the expectation expressed in Ref. [20| that 



sensitive tests of QCD in low energy vrvr scattering should use all theoretical constraints and pertinent 
low energy observables. It should, however, be emphasized that we have considered only quantities that 
do not require vital chiral inputs such as a^: work is in progress to extract sharp predictions for these 
quantities. 

After this work was completed, we have received a preprint [22] where several O(p^) coupling 
constants have been evaluated. It would be important to study the implications of this evaluation to the 
program discussed here and vice versa. 

Acknowledgements: It is a pleasure to thank H. Leutwyler for discussions and insights. BA thanks 
B. Moussallam and PB thanks J. Gasser for discussions. 



A Derivation of the Sum rules 

The three lowest partial waves can be written as: 

Re/O^ = a[] + (2a[l-5a2)^ 

+4EE^/ du'K^^^'{u,u')lmff{u') (9) 

I'=Ol'=0 "'0 

Refliu) = (2a°-5a2)^ 
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2 oo 



+4 E E ^ / du'K({{u, u')lm fl'iu') (10) 

I'=Ol'=0 

Re/o'(^.) = a2-(2a[]-5ag)^ 

+4 E E ^ / ^')Im /.^(z.') (11) 

As an example we give the derivation of the higher threshold parameters for I = 0, / = 0. 
Working in the S- and P-wave approximation, the summation over the angular momentum I runs from 
zero to one. Besides, this is equivalent to projecting eq.(^ on the relevant partial wave. The contributions 
of the three lowest partial waves to /g (j^) are explicitly written down as: 



fiv) = -P du'l ^ - + — In^ —}lmf^(u 

fit. IpP^'J 3(31/ + 2;.^ + 4) , Q{2u + u' + 2) ^ y + ^' + i ^ 

f {u) = -P du \ + In \lmf^{u), (12) 

TT Jo [ + V') VV' U' + 1 J 

/H = -P du' I \ +— In ——\\mf^{v'), 

vr Jo y6u'{l + v') 3u v' + l } 

such that 

Re /°(z.) = al + (2al - Sa^)^ + f{u) + f\u) + f\v). (13) 

The singularity of the integrand resides in the first term of f^{y). Adding and subtracting Im /q (z^') / (vr(z^+ 
z^' + 1)) does not change the integral and we may write 

f{v) = G{u) + 

I r , ,( 2 u + u' + l 2 3 + 5u' 1 \, ror 

- du'\—ln — „ ^ - „ ^ - ^ Im /o" z/tM) 

TT Jo Vov v' + 1 ^v'il + v') 3i/'(l + z^'j z^ + z^' + 1 J 

where 
Consider the difference 



1 r°° -I- 1 



f»-f» = {2al-half- + G{v)-G{Q) 

+f{^) - fiO) + fHi^) - f\0) + fii^) - f{0) (16) 
f{u) = /°(z.)-G(z.) 
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We note that the integrals not involving G in eq. ( |16| ) are free of singularities. Furthermore G(z^) — G(0) 
is also free of singularities as we show below. Consider 

G(.) - G(0) = IP /" -J>^±M±jf),r, /«(„') (17) 



which still contains the Principal Value singularity. With 

dr' 



r'2-T2 



(18) 



and t' = v'{v' + 1) and r = v{v + 1) we can write 



P / dv' , =- — — = 0. (19) 

Jo ^v'{v' + l){v> -v){v' + V+1) 

Furthermore, in the S- (and P-) wave approximation, assuming normal threshold behavior, we have: 

1 



\rrif^{v) = —^u{u+l)a\u) (20) 
and therefore we can write 

_ G(o) = -L r ^'"^r^j [^'w - ^»} (2.) 

4vr^ vz/'(z^' + l)(z^' - + I/' + 1) I- J 

which is seen to be free from the Principal Value singularity. 

The higher threshold parameters are now given by the Taylor coefficients of the difference /q (i^) — 
/o°(0)- 



B Sum rules for the higher threshold parameters 



For the higher threshold parameters considered in the text and which have not been considered in the 



literature [21 1, we find 



eg 



64 /-"^ , 
— / dv 
vr Jo 



+ 



1 



+ 



288(1 + i/)3 
1 



64i/3 288(1 + 



Im/o'(i^) 



l + 2v 



2567rz^3(l + z^)3 



(1 +2z^)x 



(22) 
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256 r 
vr Jo 



+ 



dv 



1536(1 + 
1 



512z^(l + 
Vz^(l + i^) 



2561/4 1536(1 + ^ 10247rz/(z/ + 1) 

[ l + z/ + z/^ d 0. M , l + 2i/ + 2z/2 0^^^ , 
X^^HTTWdi^" ^'^''^=° ^ .3(1 + ^)3 " + 

Similar one gets the threshold parameters for other waves: 



(1 + 2u) 
1 <P 

2{U + dzy2 



■<T\y)\u=o 



256 r 



eg 



+ 



64 



TT Jo 



2560(1 + z/)4 
2 + llzy 



1024(1 + 



rlm/o'(i/) 



2567r 5120(1 + v) 



T A( ^ V^^(l + z^) 1 + 2;/ ^ 



du 



1 



+ 



+ 



288(1 + z/)3 
1 



Im/0(i/) + 



30727rz/(l + z/) 1/ + ^^,2 



64z/(l + i^] 



Mi^ ' 576(1 + z/)3/"'^o(^^ 2567rl3'(l'+i)3^^ + ^''^'' 



dl 



|(1 + 1/ + Z/2)a2(z/) + (z/ + j.2)i_^2(^)|^^^ 

1 



2 + 5z/ 



+ 



1024z/(l + u) 



256i/4 3072(1 + 

1 + Z/ + Z/2 d 



l + 2z/ + 2z/2 1 ^2 y 

.2(i + .)3 " («) + 5(^;T^d;;2- Mi-o|^ 



\z/2(l + j,)2rfj. 

In the same way one gets the threshold parameters for / = 0, Z = 2, / = 2, Z = 2: 



hi 



dv- 



SOtt 7o '^(l + V 



dv- 



607r Jo ^(1 + ^) 



-z/Im /0(i/) + (31/ - 6)Im (i/) - 5i/Im /^(i/) 
-2i/Im /0(i/) - (3z/ - 6)Im /^(z/) - v\m f^{u) 



(23) 



(24) 



(25) 



(26) 



(27) 
(28) 
(29) 



A check on the sum rules for d^, cl and 62 is to saturate the right hand sides with the lowest order 
chiral phase shifts which should then yield the pure one-loop formulae for these quantities due to the 
perturbative unitarity of the chiral expansion. The sum rules for Cq do not converge fast enough; one 
may consider appropriate hnear combinations with say the D- wave scattering lengths in which the 
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h,2 dependence drop out and to saturate the corresponding sum rules with the lowest order chiral phase 
shifts and reproduce the one-loop formulae for these combinations. 
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44- 

if 




% 


r,2 
«0 


?,2 

Oo 


«i 


Ol 


I 


0.19 


0.237 


-0.04 


-0.074 


0.035 


0.006 


II 


0.20 


0.236 


-0.036 


-0.074 


0.035 


0.006 


III 


0.21 


0.236 


-0.035 


-0.075 


0.035 


0.006 


IV 


0.26 


0.231 


-0.021 


-0.076 


0.036 


0.006 



Table 1. List of inputs for the effective range formula 



# 


Cq 


2 
Co 


"o 


dl 


cl 




bl 


al 


I 


0.0084 


0.012 


-0.016 


-0.0045 


0.00080 


-0.00031 


-0.00028 


0.000049 


II 


0.0061 


0.013 


-0.016 


-0.0047 


0.00084 


-0.00031 


-0.00029 


0.000051 


III 


0.0041 


0.013 


-0.016 


-0.0050 


0.00087 


-0.00032 


-0.00030 


0.000053 


IV 


-0.0049 


0.015 


-0.017 


-0.0061 


0.00096 


-0.00036 


-0.00035 


0.000060 


/2 


0.0016 


0.00024 


-0.000030 


-0.0000304 


5.7- IQ-s 


0.000016 


2.2- IQ-s 


3.2- IQ-'^ 



Table 2. Higher threshold parameters computed from sum rules for the inputs of Table I. The 

contribution of /2 is also listed. 



# 


Co 


Co 


4 




Cl 


60 


bl 


al 


0-L 


0.0075 
0.0081 
0.0087 


0.015 
0.015 
0.015 


-0.0066 


-0.0036 


-0.00028 


-0.00039 


-0.00022 


0.000020 


T-L(S) 


0.0068 


0.013 


-0.017 


-0.0049 


0.00032 


-0.00031 


-0.00031 


0.000050 


T-L(G) 










0.00046 


-0.00034 


-0.00035 


0.000058 



Table 3. Table of values of quantities computed at one- and two-loop (standard and generalized) order. 
For the C/ the quantities are computed from Zi and I2 listed in Table 1 of Ref. [13]. The two-loop 

quantities are from [pi] and [IT] 
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